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Abstract. We revisit the determination of as{m^) using a fit to inclusive r hadronic spectral 
moments in light of (1) the recent calculation of the fourth-order perturbative coefficient K4 
in the expansion of the Adler function, (2) new precision measurements from BABAR of 
e^e~ annihilation cross sections, which decrease the uncertainty in the separation of vector 
and axial-vector spectral functions, and (3) improved results from BABAR and Belle on r 
branching fractions involving kaons. We estimate that the fourth-order perturbative predic- 
tion reduces the theoretical uncertainty, introduced by the truncation of the series, by 20% 
with respect to earlier determinations. We discuss to some detail the perturbative prediction 
of two different methods: fixed-order perturbation theory (FOPT) and contour-improved 
perturbative theory (CIPT). The corresponding theoretical uncertainties are studied at the 
T and Z mass scales. The CIPT method is found to be more stable with respect to the miss- 
ing higher order contributions and to renormalisation scale variations. It is also shown that 
FOPT suffers from convergence problems along the complex integration contour. Nonpertur- 
bative contributions extracted from the most inclusive fit are small, in agreement with earlier 
determinations. Systematic efi'ects from quark-hadron duality violation are estimated with 
simple models and found to be within the quoted systematic errors. The fit based on CIPT 
gives asim^) — 0.344 ± 0.005 ± 0.007, where the first error is experimental and the second 
theoretical. After evolution to Mz we obtain as(M|) = 0.1212 ± 0.0005 ± 0.0008 ± 0.0005, 
where the errors are respectively experimental, theoretical and due to the evolution. The 
result is in agreement with the corresponding N'^LO value derived from essentially the Z 
width in the global electroweak fit. The as{Mz) determination from t decays is the most 
precise one to date. 
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1 Introduction 

The relatively large mass of the r lepton, its leptonic nature and its decay through weak interaction 
promotes it to a particular status for probing the Standard Model (see [1] for a detailed review, and 
references therein). In particular, spectral functions determined from the invariant mass distributions of 
hadronic r decays are fundamental quantities describing the production of hadrons from the non-trivial 
vacuum of strong interactions. They embed similar information to the one determined from cross sections 
of e^e~ annihilation to hadrons: both kinds of spectral functions are especially useful at low energies 
where perturbative QCD fails to locally describe the data, and where the theoretical understanding of the 
strong interactions remains at a qualitative level. Due to these limitations on the theoretical side, spectral 
functions play a crucial role in calculations of hadronic vacuum polarisation contributions to observables 
such as the effective electromagnetic coupling at the Z mass, and the muon anomalous magnetic moment. 

Inclusive hadronic quantities, obtained after integrating over the spectral functions (or directly via the 
measurement of hadronic or leptonic r branching fractions) , have been found to be dominated by pertur- 
bative contributions at energies above ~1 GeV. They can be exploited to precisely determine the strong 



2 



coupling constant at the r-mass scale, as(m^) [2, 3, 4, 5]. More recently, this determination was re- 
assessed [1] in the light of the existing data on t decays and e''"e~ annihilation. 

In the present paper, we update the determination of as(m^) from hadronic r decays, motivated by 
progress performed in two different areas: on the theoretical side, the perturbative expression of the rel- 
evant correlator has been computed up to fourth order [6], and on the experimental side, new precision 
measurements from BABAR of t branching fractions involving kaons [7] decrease the uncertainty in the 
separation of vector and axial-vector spectral functions. We utilise this opportunity to analyse several 
features of the theoretical frameworks commonly used to determine asirn^) in more detail. This concerns 
the treatments of the perturbative series, the convergence of the expansions, and the impact of nonper- 
turbative effects. 

In Sec. 2 we describe recent experimental improvement on the measurements of KKt: decays, the spec- 
tral functions and the t branching fractions. This is followed in Sec. 3 by a summary of the various 
theoretical prescriptions used to extract as{iTi^) from a fit to data, and a discussion of their advantages 
and shortcomings. We also analyse the role played by nonperturbative contributions in this determina- 
tion. In Sec. 4 we exploit the normalisation and shape of the spectral functions to constrain the relevant 
nonperturbative contributions and to provide an improved determination of as(m^). 



2 Tau Hadronic Spectral Functions 

For vector (axial- vector) hadronic r decay channels V^v-r (A^Vr), the nonstrange vector (axial- vector) 
spectral function vi (ai, oq), where the subscript refers to the spin J of the hadronic system, is derived 
from the invariant mass-squared distribution (\/Ny/A){dNy/A/ ds) for a given hadronic mass y/s, divided 
by the appropriate kinematic factor, and normalised to the hadronic branching fraction 
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For ao(s), the same expression holds if the term (1 + 2s/m^) is removed. Here S'ew = 1.0198 ± 0.0006 is a 
short-distance electroweak correction [8, 9], By- /a-v^ i^e) denotes the inclusive r /A~ {j)^^ (r 

e~Fei^r) branching fraction (throughout this letter, final state photon radiation is accounted for in the t 
branching fractions) . We use universality in the leptonic weak charged currents and the measurements of 
Be, Bf, and the r hfetime, to obtain the improved branching fraction Be = Bf' = (17.818±0.032)% [1]. We 
also use = (1776.90 ± 0.20) MeV [10] and \Vud\ = 0.97418 ± 0.00019 [11] (assuming CKM unitarity). 
Integration of the spectral function over the r phase space leads to the inclusive r hadronic width, 
expressed through the ratio 
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By unitarity and analyticity the spectral functions are connected to the imaginary part of the two-point 
correlation function, n^^jj{q), for time-like momenta-squared > 0, 

n^;',j{q) ^ ^ J d^xe^^^{0\TiUt^ixM^iOnO) = i-g^-'q'+q^qn n\%{q^) + q^^ q'^ n\%{q^) , (3) 

where U — A,V denotes the nature of the relevant currents, either vector {Ul^j — V^^ — q.p^qi) or 
axial-vector {U^, — A^- = qj^f^^^qi) charged colour-singlet quark currents. By Lorentz decomposition, 
the correlation functions can be split into their J = 1 and J = parts. 

In the complex s — q^ plane, the polarisation functions II^^jj{s) are expected to exhibit a very simple 
analytic structure, the only non-analytic features being along the real axis: a branch cut for all polarisation 
functions, and a pole at the pion (kaon) mass for oq. The imaginary part of the polarisation functions 
on the branch cut is linked to the spectral functions defined in Eq. (1), for nonstrange (strange) quark 
currents ^ 

Im-^4dS,yM(«) = 2^^^i/«i,o(s) , (4) 
which provide the basis for comparing a theoretical description of strong interaction with hadronic data. 
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Experimentally, the total hadronic observable Rt, 

Rr Rr.V + Rr,A + Rr,S , (5) 

where Rt,s denotes the hadronic width to final states with net strangeness, is obtained from the measured 
leptonic branching ratios, 

^ 1 _ (6) 

Be 

2.1 New Input to the Vector/ Axial- Vector Separation 

The separation of vector and axial-vector components is straightforward in the case of hadronic final 
states with only pions using G-parity, provided that isospin symmetry holds. An even number of pions 
has G — 1 corresponding to vector states, while an odd number of pions has G — —1, which tags axial- 
vector states. Modes with a KK pair are not in general eigenstates of G-parity and contribute to both V 
and A channels. While the decay to K~K^ is pure vector, additional information is required to separate 
the KKiT and the rarer KKtttt modes. For the latter channel an axial-vector fraction of 0.5 ± 0.5 is 
used [1]. 

Until recently, there was some confusion on this issue for the KKtt modes: 

1. In the ALEPH analysis of r decay modes with kaons [12], an estimate of the vector contribution 
was obtained using the e+e^ annihilation data from DM1 [13] and DM2 [14] in the KKtt channel, 
extracted in the 1=1 state. This contribution was found to be small, and, using the conserved vector 
current (CVC), a branching fraction of Bcvc{t Vt{KKi:)v) = (0.26 ± 0.39) • 10"'^, was found, 
corresponding to an axial fraction of fA,cvc{KKiT) = 0.94to o8- 

2. The ALEPH CVC result was corroborated by a partial-wave and lineshape analysis of the ai resonance 
from T decays in the Vt-'k~'2,'k^ mode performed by CLEO [15]. The effect of the K*K decay mode of the 
fli was seen through unitarity and a branching fraction of B{ai K*K) = (3.3 ± 0.5)% was derived. 
With the known i/r^i branching fraction, this value more than saturates the total branching 
fraction available for the KKtt channel, yielding an axial fraction of fA,ai {KKtt) = 1.30 ± 0.24. 

3. Another piece of information, also contributed by CLEO [16], but conflicting with the two previ- 
ous results, is based on a partial- wave analysis in the K~ K^tt~ channel using two-body resonance 
production and including many possible contributing channels. A much smaller axial fraction of 

KKTri^^^) — ^-^^ ^ 0.10 was found here. 

Since the three determinations are inconsistent, the value — 0.75 ± 0.25 has been used previously 
to account for the discrepancy [1]. This led to a systematic uncertainty in the V, A spectral functions that 
competed with the purely experimental uncertainties. 

Precise cross section measurements for e+e^ annihilation to K^ K^tt^ and to K^K^tt^ have been recently 
published by the BABAR Collaboration [7], using the method of radiative return. In the mass range of 
interest for r physics they show strong dominance of K*{890)K dynamics and a fit of the Dalitz plot 
yields a clean separation of the / = 0, 1 contributions. Assuming CVC, the mass distribution of the vector 
final state in the decays r —>■ Vt-KKtt can be obtained. The result is shown in Fig. 1 and compared with 
the full T spectrum from ALEPH [12] summing up the contributions from the K~ K^tt~ , K^K^tt~, and 
K~K'^tt'^ modes. The BABAR results reveal a small vector component. After integration, one obtains 

/a,cvc {KKtt) = 0.833 ± 0.024, (7) 

which is about 1.3(t lower than the ALEPH determination using the same method (but with much poorer 
e~^e~ input data) and 2.7f7 higher than the CLEO partial- wave-analysis result. The new determination 
has a precision that exceeds the previously used value by an order of magnitude, thus effectively reducing 
the uncertainties in the vector and axial-vector spectral functions to the experimental errors only. 

One notices from Fig. 1 that the axial fraction varies versus the KKtt mass, with lower masses being 
further axial-enhanced. The observed axial- vector dominance is at variance with several estimates such as 
/a ^ 0.10 [17], 0.37 [18], obtained within the Resonance Chiral Theory, which attempts at incorporating 
massive vector and axial resonances decaying into light mesons into a framework inspired by chiral and 
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Fig. 1. The mass-squared distribution for r Vt-KKtv decay modes from ALEPH and the predictions for the 
vector component obtained by CVC using DM1, DM2 and BABAR e+e" data. 



large- A^c arguments. On the other hand, this axial-vector dominance is closer to the prediction fj^ ~ 0.71, 
based on a model combining axial- vector and vector resonances of finite widths with a leading-order chiral 
Lagrangian [19]. 

In deriving Eq. (7) care was taken to include a small contribution from the (f>TT final state, observed 
by BABAR in the same analysis [7]. Since BABAR also published a ^ ly-r'^'n'^ branching fraction 
measurement [20] , it is possible to perform a test of CVC in this channel with 



Bcvcir Vr(j)TT^) = (3.8 ± 0.9 ± 0.2) • 10 
Br{T ^ l^r^TT^) = (3.42 ± 0.55 ± 0.25) ■ 
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for which we find agreement within the quoted statistical and systematic errors. For comparison the 
dominant CVC r -> i^rK*{890)K branching fraction is (7.3 ± 0.6 ± 0.4) • lO""*. 



2.2 Update on the Branching Fraction for Strange Decays 

New measurements of r strange decays have been published since our last compilation [1]. This is the case 
for the hadronic channels Kn'^ [21], Ks7r~ [22], and K~TT~^Tr~ [7]. Also using the more precise estimate 
from universality for the channel [1], the updated value of i?r,s becomes 

s = 0.1615 ±0.0040, (10) 

replacing the previous value of 0.1666 ± 0.0048 [1]. 

Using the new fAiKKir) value (7), the updated hadronic widths Rr,v/A from ALEPH, slightly renor- 
malised so that their sum agrees with the new average for Rt,v+a obtained from (6) and (10) read 

R^ y = 1.783 ± 0.011 ± 0.002 , (11) 

Rr,A = 1.695 ± 0.011 ± 0.002 , (12) 

ii'r,y+A = 3.479 ±0.011 , (13) 

i?r,y-A = 0.087 ±0.018 ±0.003, (14) 

where the first errors are experimental and the second due to the V/A separation, now dominated by the 
KKirn channel. 

The ALEPH spectral functions are updated accordingly and shown in Fig. 2 for respectively vector, 
axial-vector, V + A and V — A. 
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Fig. 2. Vector (V), axial- vector (A), V + A and V — A t hadronic spectral functions measured by ALEPH, and 
updated using the new V, A separation in the KK-k channels discussed in the text. The shaded areas indicate the 
main contributing exclusive t decay channels. The curves show the predictions from the parton model (dotted) 
and from massless perturbative QCD using as (Mz) = 0.120 (solid). 



3 Theoretical Prediction of Rt 



Tests of QCD and the precise measurement of the strong coupling constant as at the r mass scale [2, 3, 
4, 5], carried out first by the ALEPH [23] and CLEO [24] collaborations, have triggered many theoretical 
developments. They concern primarily the perturbative expansion for which different optimised rules have 
been suggested. Among these are contour-improved (resummed) fixed-order perturbation theory [26, 27], 
effective charge and minimal sensitivity schemes [28, 29, 30, 31, 32], the large-/3o expansion [33, 34, 35], 
as well as combinations of these approaches. Their main differences lie in how they deal with the fact 
that the perturbative series is truncated at an order where the missing part is not expected to be small. 
While a review and discussion of the various approaches can be found in [1] , we only recall some of their 
salient features in the following. 

With the publication of the full vector and axial-vector spectral functions by ALEPH [36, 37] and 
OPAL [25] it became possible to directly study the nonperturbative properties of QCD through V — A 
sum rules and through fits to spectral moments computed from weighted integrals over the spectral func- 
tions (we refer again to the discussions in [1]). Inclusive observables like Rr can be accurately predicted 
in terms of as(TO^) using perturbative QCD, and including small nonperturbative contributions within 
the framework of the Operator Product Expansion (OPE) [38]. 
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3.1 Operator Product Expansion 

According to Eq. (4) , the absorptive (imaginary) parts of the vector and axial- vector two-point correlation 
functions n^^Jy^^(s), with the spin J of the hadronic system, are proportional to the r hadronic spectral 
functions with corresponding quantum numbers. The nonstrange ratio R^.v+A can be written as an 
integral of these spectral functions over the invariant mass-squared s of the final state hadrons [4] 



so 

f ds f s 

Rr,V+A{so) = 127rS'Ew|KdP / — — 

J So \ So 





2 

1 + 2— ] Imij(i) (s + is) + ImTjf") (s + ie) 
•So 



(15) 



where U^'^'^ can be decomposed as U^'^'^ = H^^y + h'^^a- work in the chiral limit^ to study the 
perturbative contribution, so that the lower integration limit is zero because of the pion pole at zero 
mass. The correlation function 77*^''^ is analytic in the complex s plane everywhere except on the positive 
real axis where singularities exist. Hence by Cauchy's theorem, the imaginary part of U^"^^ is proportional 
to the discontinuity across the positive real axis, and the integral (15) can be replaced by a contour integral 
over n{s) running counter-clockwise around the circle from s = sq + ie to s = So — i£- 

The energy scale sq = is large enough that contributions from nonperturbative effects are expected 
to be subdominant and the use of the Operator Product Expansion is appropriate. The latter is expected 
to yield relevant results in the deep Euclidean region where s is large and negative, whereas the extension 
to other regions in the complex plane is questionable. Fortunately, in the case of the kinematic factor 
(1 — s/so)^ suppresses the contribution from the region near the positive real axis where n^'^\s) has a 
branch cut and the validity of the OPE is doubtful due to large quark-hadron duality violations [39, 40]. 

The OPE of the vector and axial- vector ratio Rr,v/A can be written as 



Rr,V/A - ^5ew|K.P ( 1 + S^"^ + 5W + ^-.y/i + ECVm ) ' (16) 
^ D=4,6,... ^ 

with the massless universal^ perturbative contribution S^^\ the residual non-logarithmic electroweak 
correction S^y^ — 0.0010 [43] (c/. the discussion on radiative corrections in [1]), and the dimension D ^ 2 

perturbative contribution from massive quarks. The term (5*^^^ denotes the OPE contributions of 

mass dimension D [5] 

^ud'v/A = ^V/A (so , /i) , (17) 

dimO=_D ■^0 

where ~ li^ud^v + ^ud\)- practice, the OPE provides a separation between short and long 

distances by following the flow of a large incoming momentum. The scale parameter /x separates the long- 
distance nonperturbative effects, absorbed into the vacuum expectation value of the operators (©^(m)): 
from the short-distance effects that are included in the coefficients Cv/a{s, Ij), which become Cy^^(so, 
after performing the integration (15). The vacuum expectation values {OoifJ-)) encode information on 
the nonperturbative features of QCD vacuum and its effects on the propagation of quarks: they cannot 
be computed from first principles and have to be extracted from data. The short-distance coefficients 
Cy/ji{s, fi) can be determined within perturbative QCD. 



3.2 Perturbative Contribution to Fourth Order in as 

Rt is a doubly inclusive observable since it is the result of an integration over all hadronic final states at a 
given invariant mass and further over all masses between ttItt and rrir- The scale rrir lies in a compromise 
region where as(m^) is large enough so that Rr is sensitive to its value, yet still small enough so that 
the perturbative expansion converges safely and nonperturbative power terms are small. The prediction 
for Rr is thus found to be dominated by the lowest-dimension term in Eq. (17), i.e., the term obtained 

^ Vector and axial-vector currents are conserved in the chiral limit, so that sUy^ = siT^'' = 0. 
In the chiral limit of vanishing quark masses the contributions from vector and axial-vector currents coincide 
to any given order of perturbation theory and the results are flavour independent. 
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from a perturbative computation of the correlator 77. 

For the evaluation of the perturbative series, it is convenient to introduce the analytic Adler function [45] 
D(s) = —s-dn(s)/ds, which avoids extra subtractions that are unrelated to QCD dynamics. The function 
D(s) calculated in perturbative QCD within the MS renormalisation scheme is a function of as and 
depends on the renormalisation scale /x, occurring through ln(/i^/s). Since D{s) is connected to a physical 
quantity, the spectral function Im7T(s), it cannot depend on the choice of the renormalisation scale fi. 
This is achieved through the cancellation of the /i-dependence of as and of the explicit occurrences of 
H in D. Nevertheless, in the realistic case of a series truncated at a given order in as our knowledge 
of the renormalisation scale dependence is imperfect, i.e., D depends on ^, thus inducing a systematic 
uncertainty. 

To introduce the Adler function in Eq. (15), one uses partial integration, giving 

1 + ^(0) = -2^i ^ ^u;(s)7?(s) , (18) 

\s\ = so 

where w{s) — 1 — 2s/ sq + 2{s/so)^ — (s/sq)*. The perturbative expansion of D{s) reads 

D{s) = ^f2Kn{0<{-^s), (19) 

n=0 

with Us = as/TT, and where the dimensionless factor ^ parametrises the renormalisation scale ambiguity. 
While the coefficients 7fo,i = Kq^i = 1 are universal (we use the notation Kn = Kn{£, — 1) in the 
following), the K„>2 depend on the renormalisation scheme and scale used. Powerful computational 
techniques have recently allowed to determine K4. The authors of [6] exploited the dependence of the 
four-loop master integrals (used to express all relevant four-loop integrals with massless propagators) on 
the space-time dimension to compute the integrals to the required accuracy. For Uf — 3 quark flavours 
and ^ = 1 one has^ K2 ~ 1.640, K3 ~ 6.371 and K4 ~ 49.08 [46, 47, 48, 49, 50, 6]. The fuU expressions 
for the functions Kn (C) for arbitrary ^ up to order n — 5 can be found in [1] . 

With the series (19), inserted into the r.h.s. of Eq. (18), one obtains the perturbative expansion 

oc 

<5(") -^7^„(0A(")(a,), (20) 

n=l 

with the functions [26] 

TT 

^^"^as) <f -w{s)al^i-^s) = ^ / dipwi-SoeV^i^soe^n ■ (21) 

zni J s 271 J 



Similarly, the Adler function also serves to obtain the perturbative expansion of the inclusive e+e 
annihilation cross section ratio 

(T(e+e^ ^ hadrons (7)) -V^r^a / j / ^i^') 



Re+e- (s) = , ^ _ = -dm ^ Qj f ds' ■ (22) 



cr(e+e fi^ii ) 



Evaluating the contour integral in fixed-order perturbation theory (c/. Sec. 3.2.1) with nf = 5 active 
quark flavours, and inserting all known coefficients, gives'* 



Rf+^-{s) = 3^^gJ [l + a,(s) + 1.4092 a2(s) - 12.7673 a3(s) - 79.9795 a^(s 

+ (7^5 + 79.7306) al{s) + {Kq + 2202.78) al{s) + ...]. (23) 



^ The numerical expressions for an arbitrary number of quark flavours (n/) in the MS renormalisation scheme 
for C = 1 are: Kq ^ 1, Ki ^ 1, K2 ~ 1.9857 - 0.1153n/, K3 ~ 18.2428 - 4.2158n/ + 0.0862ny, and K4 ~ 
135.7916 - 34.4402 n/ + 1.8753 n} - 0.0101 n}. 

* The explicit formula reads: 
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3.2.1 Fixed-Order and Contour-Improved Perturbation Theory 

The standard perturbative method to compute the contour integral consists of expanding all the quantities 
up to a given power of as(so)- The starting point is the solution of the renormalisation group equation 
(RGE) for as{s), which is expanded in a Taylor series of 77 = ln(s/so) around the reference scale sq [1] 



3.3 \ 4 



+ + \pW + 3/30/327?' - ^PlPirf + PW^ a' (24) 

+ - /3477 + \piP2V^ + ^/3o/33^' - - 6/3o/327?' + ^^HlPiV^ - /3oS')«s + 0{v';al) . 

Here the series has been reordered in powers of = as(so) and we use the RGE /3-function^ as defined 
in [51]. 

Computing the contour integral (21), and ordering the contributions according to their powers in as, 
leads to the familiar expression for fixed-order perturbation theory (FOPT) [26] 

00 

^^°^ = E[^«(^)+5n(0]«"(^so), (25) 

n=l 

where the (?„ are functions of Km<n and Pm<n-i, and of elementary integrals with logarithms of power 
m < n in the integrand. Setting = 1 and replacing all known (3i coefhcients by their numerical values 
for ri/ = 3 gives [1, 52] 

,5(0) = as{so) + {K2 + 3.5625) al{so) + {K3 + 19.995) al{so) + {K^ + 78.003) a^(so) (26) 
+ {K5 + 307.787) a^(so) + (Ke + 17.813 A'5 -I- 1.5833 - 5848.19) a%so) , 

where for the purpose of later studies we have kept terms up to sixth order. 

The FOPT series is truncated at a given order despite the fact that parts of the higher coefficients 
5n>4(0 are known and could be resummed: these are the higher order terms of the as(s) expansion that 
are functions of /3„<3 and KnKi only. Moreover, at each integration step, the expansion (24) with respect 
to the physical value as(so) is used to predict as{s) on the entire |s| — sq contour. This might not always 
be justified, and leads to systematic errors as discussed in Sec. 3.2.3. 

A more accurate approach to the solution of the contour integral (21) is to perform a direct numerical 
evaluation by step-wise integration. At each integration step, it takes as input for the running as(s) the 
solution of the RGE to four loops, computed using the value from the previous step [27, 26]. It implicitly 
provides a partial resummation of the (known) higher order logarithmic contributions, and does not 
require the validity of the as{s) Taylor series for large absolute values of the expansion parameter rj. This 
numerical solution of Eq. (20) is referred to as contour-improved perturbation theory (CIPT). 



3.2.2 Alternative Perturbative Expansions 

Inspired by the pioneering work in [28, 29, 30, 31, 32] the effective charge approach to the perturba- 
tive prediction of Rt (ECPT) has triggered many studies [53, 54, 55]. The advocated advantage of this 
technique is that the perturbative prediction of the effective charge is renormalisation scheme and scale 



^ The full expressions for an arbitrary number of quark flavours (n/) are: /3o = j (11 — f'l-/), Pi = 

1 /ino 38^ \ _ 1 {2S57 5033 „ ,325 2\ _ 1 [ 149753 , ocfi/l /- / 1078361 , 6508 a \ „ , 

16 (102- —Uf), 1^2 - {— jg-n/ -I- —Uf), and Ih - [— g— +'^564C3 - ( lea + — Caj + 

^52065 _^ 6472 ^2 ^ l^n)] , where the ^={3,4,5} = {1.2020569, vr^/gO, 1.0369278} are the Riemann C-functions. 
The /3„>4 are unknown. 
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invariant since it is a physical observable. The effective r charge is defined by ar = S'^^\ The ECPT 
scheme has been used in the past to estimate the unknown higher-order perturbative coefRcient K^, 
by exploiting the mediocre convergence of the series (because arim^) ~ f .8 • as(m^)). As pointed out 
in [6], these estimates missed the actual value of K4 by approximately a factor of two. One reason for 
this disagreement may come from the fact that these methods neglected the contributions from the next 
higher and also unknown orders. Owing to the insufficient convergence, the uncertainty on the coefficient 
estimate introduced by this neglect is significant and exceeds the errors quoted [1] . 

For completeness we also mention the large-(3o expansion, which is an approximation to the full FOPT 
result assuming the dominance of the [Poas^—s)]" term. It is thus possible to derive estimates for the 
FOPT coefficients of a given perturbative series at all orders by neglecting higher order terms in the 
/3-function. The large-/3o expansion corresponds to inserting chains of fermion loops into the giuon prop- 
agators and to determining the impact on the quark- antiquark vacuum polarisation. The procedure 
provides hence a naive non-abelianisation of the theory, because the lowest-order radiative corrections 
do not include gluon self-coupling. As an illustration, the Rr FOPT series (25) can be expanded as 
5(o)(s) = Us J2n=o '^s {dnPo + Sn) , where dnPo + (S„ = Kn+i + gn+i (setting ^ = 1). The coefHcients d„ 
are computed in terms of fermion bubble diagrams [56] , where they are identified with their leading-n / 

pieces Sn^ in the expression dn = d^n^n'j: + . . . + d'^^\ Neglecting the corrections (5„, the above series leads to 

the large-/?o expansion of S''^\ The first elements of the series are [57]: do = 1, diPo = 5.119, (i2/9o = 28.78, 
^3/33 = 156.6, d4/3* = 900.9, ds/?^ = 4867. They compare reasonably well with the FOPT terms (26) where 
these are known, in particular the large size of the fourth-order term has been anticipated {K4 ~ 79). 
However, it turns out that the estimated coefficients of the Adler series itself (before integration on the 
contour) do not compare well with the exact solutions, which emphasises the uncontrolled theoretical 
uncertainties associated with this method [1] . 

3.2.3 Comparing Perturbative Methods 

This section updates and completes the discussion given in Sees. 3 and 8 of [1], including here the known 
value of the fourth-order perturbative coefRcient in the Adler function, K4 [6]. We perform a numerical 
study of the FOPT and CIPT approaches to expose the differences between these two methods. Both use 
the Taylor series (24), and they assume that one can perform an analytic continuation of the solution of 
the RGE for complex values of s,^ namely along the circular contour of integration in Eq. (21). One should 
thus make sure that the series is used only inside the domain of good convergence. As one approaches the 
limit of this domain, the error induced by the finite Taylor series increases. For CIPT the convergence is 
guaranteed because the integration proceeds along infinitesimal steps such that |r/| ^ 1 everywhere. The 
situation is more complicated for FOPT as the absolute value of 77 in Eq. (24) approaches tt close to the 
branch cut. 

The tests carried out here use the expansion (24) to sixth order in as (so) (hence fifth order in 77 = 
ln(s/so)) — if not stated otherwise, with estimates for K^ q and (i^ assuming a geometric growth of the 
corresponding series {i.e., i^5(6) = -f^4(5)(J^4(5)/^3(4)) and /34 = PsiPs/ 132)), and setting all coefficients 
at higher-orders than these to zero. 



Taylor Series 

To check the stability of the results obtained with FOPT, we consider a variant (denoted FOPT^^) where 
all known or estimated terms of order 77"-^ are kept (i.e., including the known expressions with powers 
a"^^(so) and beyond), which should reduce the error associated with the use of the Taylor expansion in 
FOPT. Figure 3 shows the evolution of the real part of as(s) along the integration circle as found for 
CIPT, FOPT and FOPT++. As expected, the values for CIPT and F0PT(++^ agree in the region around 
= (the fix-point of the expansion in FOPT'-"''^^), but significant discrepancies occur elsewhere. For 
FOPT'''"^ we find large values for Re(Q;s) close to the branch cut. Estimating the convergence speed of 
the 77 series (24) reveals that it is slower for FOPT"*""*", where larger powers of are kept, than for FOPT, 
for which the series is truncated at . Including higher orders ?7">5 in FOPT"'"'' we find that these terms 

® One of the first limits of this hypothesis shows up in the discontinuity of the imaginary part of as ai (f> = ±7r, 
which is due to the cut of the logarithm in the complex plane. 
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Fig. 3. Real part of as{s) computed along the \s\ = so contour for ^ = 1, using respectively FOPT"*"^ (solid line, 
see text), FOPT and FOPT+ (dashed, see text) and CIPT (dashed-dotted) . 
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Fig. 4. Real part of (47r^D(s) — 1) (left) and of the integrand in Eqs. (21) and (20) (right), computed along the 
integration contour for ^ = 1, using respectively FOPT++ (solid line), FOPT (dashed), CIPT (dashed-dotted) 
and FOPT"*" (dotted, not shown on the right hand plot because it is almost indistinguishable from CIPT). 



dominate the value of Re (as) near the branch cut, leading to large deviations from the correct evolution, 
which rise with the order n. On the contrary, for CIPT performed with infinitesimal integration steps, 
the full five-loop RGE solution is equivalent to Eq. (24), i.e., CIPT = CIPT++.'^ 

Although the values of as differ significantly on half of the integration domain, the standard FOPT and 
CIPT methods give similar results for the integral. This is because the integration kernel (18) vanishes 
for s — —So (0 = iTi")! suppressing the contributions to the integral coming from the region near the 
branch cut.^ The main difference between the two results stems from the regions 4> ~ ±2.1 and cj) « ±1 
(c/. left-hand plot of Fig. 4). In the region \(f>\ < 1, the values of as{s) estimated by the two methods are 



^ To understand this feature, one can compare the errors induced by the Taylor approximation for the FOPT 
and CIPT numerical procedures along the circular contour. To compute the contour integral, N ^ 1 equidistant 
integration points along the contour are added. At the j"^ point, the error on the value of as is given directly by 
Eq. (24) for FOPT, whereas one can easily show that it is reduced by the factor j/N"'^'^ for CIPT, where n = 5 
is the expansion order in rj. Therefore, the error on the contour integral coming from the determination of as is 
suppressed by 1 /N" in the case of CIPT compared to FOPT. 

* In addition, a significant cancellation takes place in this region; for FOPT, the contribution of the contour 
integral vanishes on the intervals [— vr; —1.73] and [1.73; n], whereas for CIPT a vanishing contribution comes from 
[-7r;-1.57] and [1.57; vr]. 
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Table 1. Massless perturbative contribution S^"^ in Rt using FOPT, CIPT and the large-/3o expansion, respec- 
tively, and computed for as(m^) = 0.34. The unknown higher-order K^^e and P4 coefficients are estimated by 



assuming a geometric growth (see text), while the remaining 
correspond to the indicated error ranges. 


ones are 


set to zero. 


The quoted uncertainties 5 


Pert. Method 


n = 1 


n = 2 


n = 3 


71 = 4 


(n = 5) 


(n = 6) 


Z^n = l 


y^5 


Z^n = l 


FOPT = 1) 


0.1082 


0.0609 


0.0334 


0.0174 


0.0101 


0.0067 


0.2200 


0.2302 


0.2369 


SiPi ± 100%) 

















±0.0006 








±0.0006 


5{K5 ± 100%) 














±0.0056 


±0.0108 





±0.0056 


±0.0164 


5{K6 ± 100%) 
± 0.63) 

















±0.0047 




-1-0.0317 
-0.0151 




+0.0209 
-0.0119 


±0.0047 

+0.0152 
-0.0095 



CIPT = 1) 


0.1476 


0.0295 


0.0121 


0.0085 


0.0049 


0.0020 


0.1977 


0.2027 


0.2047 


S{f34 ± 100%) 


TO.0003 


TO.0001 


to.oooi 


TO.0001 


TO.0001 


TO.0001 


TO.0006 


TO.0007 


T0.0008 


S{K5 ± 100%) 














±0.0049 








±0.0049 


±0.0049 


5{Ke ± 100%) 

















±0.0020 


O 


O 


±0.0020 


(5(5 ± 0.63) 














+0.0032 
-0.0051 


+0.0005 
-0.0044 


+0.0001 

-0.0079 



Large-/3o expansion 0.1082 0.0600 0.0364 0.0215 0.0134 0.0078 0.2261 0.2395 0.2473 



close, and the difference between the two integrands can be ascribed to the truncation at the sixth order 
in as{so) for the integrand of FOPT. 



Fixed-order Truncation 

In addition to employing a Taylor series in a region with questionable convergence properties, FOPT 
truncates the full expression of the contour integral in Eq. (25). To disentangle the impact of these two 
approximations, we have tested another variant of FOPT (denoted FOPT"*"), where Eq. (24) is used as is, 
but without truncating the Adler function (or equivalently J^"^) at the sixth order in ag[sQ). This method 
leads to a similar integrand as in CIPT, with however the usual difference in the evolution. The left-hand 
plot of Fig. 4 shows the evolution of the real part of (47r^Z?(s) — 1) along the contour for all methods. 
FOPT+ and CIPT differ close to the branch cut as a consequence of the deficient Taylor approximation, 
with however little difference in the integration result [1] due to the suppression by the integration kernel. 
The FOPT^^ approach without truncating the Adler function leads to a J^"^ that lies between CIPT 
and FOPT, with however unstable numerical dependence on the largest power in 77 kept in the Taylor 
series. 



Numerical Comparisons 

Table 1 summarises the contributions of the orders n < 6 in PT to (5'"' for FOPT, CIPT and the large-/3o 
expansion,^ using as benchmark value Q!s(m^) — 0.34, and ^ = 1. For systematic studies we vary ^ in the 
range ^ • = ± 2GeV^ and the maximum observed deviations with respect to ^ = 1 are reported 
in the corresponding lines of Table 1. We assume a geometric growth of the perturbative terms for all 
unknown PT and RGE coefficients, with 100% uncertainty assigned to each of them for the purpose of 
illustration. We recall that the n-th contributions to the FOPT and CIPT series should be compared 
with care. Whereas the FOPT contributions can be directly obtained from Eq. (25), the entanglement 
of the different perturbative orders generated by CIPT prevents us from separating the contributions in 
powers of as{so). Instead, the columns given for CIPT in Table 1 correspond to the terms in Eq. (20). If 
the two methods were equally well suited for the integration, their column sums should converge to the 
same value. 

^ We do not include ECPT into the present study, because — as concluded in [1] — the convergence of the 
perturbative series is insufficient for a precision determination of 03(01^). 
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Fig. 5. Scale dependence of j'''' in Rt computed at the third to the estimated sixth order with FOPT (left) and 
CIPT (right). 
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Fig. 6. Scale dependence of 5^"^ in R _l _ (Aff) computed at the third to the estimated sixth order with FOPT 
(left) and CIPT (right). 



The variations of (J^*^-* with the scale parameter ^ are strongly non- linear (c/. the asymmetric errors in 
Table 1 and the functional forms plotted for FOPT (left) and CIPT (right) in Fig.' 5). CIPT exhibits 
significantly less renormalisation scale dependence than FOPT at order n = 4, while the interpretation 
of the subsequent orders strongly depends on the values used for the unknown coefficients Kn>5- 



Conclusions 



The CIPT series is found to be better behaved than FOPT and is therefore to be preferred for the 
numerical analysis of the r hadronic width. This preference is also supported by the analysis of the 
integrand in the previous section, suggesting a pathological behaviour of FOPT for Os near the branch 
cut. Our coarse extrapolation of the higher-order coefficients could indicate that minimal sensitivity is 
reached at n '--^ 5 for FOPT, while the series further converges for CIPT. The uncertainties due to and 
Kq are smaller for CIPT whereas the one due to the unknown value of (3^ is similar in both approaches. 
The difference in the result observed when using a Taylor expansion and truncating the perturbative 
series after integrating along the contour (FOPT) with the exact result at given order (CIPT) exemplifies 
the incompleteness of the perturbative series. The situation is even worse since, not only large known 
contributions are neglected in FOPT, but the series is also used in a domain where its convergence is 
not guaranteed: taking the difference between CIPT and FOPT as an estimate of the related systematic 
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Table 2. Massless perturbative contributions to i'''' in i?^+^_(Aff) using FOPT and CIPT, respectively, and 
computed for as(Af|) = 0.12. The unknown higher-order K^^ and (34 coefficients are estimated by assuming a 
geometric growth, while the others are set to 0. The quoted uncertainties S stem from the indicated range of 
values for the unknown parameters and from the renormalisation scale. 

Pert. Method n=l n = 2 n = 3 n = 4 {n = 5) {n = 6) ELi ELi ELi 

FOPT = 1) 0.038197 0.002056 -0.000712 -0.000170 -0.000004 0.000012 0.039372 0.039368 0.039380 
5(/34±100%) 000000000 
S(K5±100%) tW-^ TlO"'^ TlO"^ 

S{Ke ± 100%) ±5 ■ lO"** ±5 • 10"'' 

5(C ± 0.63) - - - - - - 10-« tit ■ 10-" tli ■ W-"" 

CIPT (5 = 1) 0.037462 0.001941 -0.000034 0.000016 -0.000008 0.000003 0.039385 0.039378 0.039381 

5{f34 ± 100%) < lO"** «0 fsO ^0 ^0 ^0< 10^^ < 10"'' < 10"^ 

5{K5 ± 100%) =f8 ■ 10"'' =f8 • lO"'' t8 • 10"" 

S{Ke ± 100%) ±3 • lO"'' ±3 • 10"'' 

<5(^ ± 0.63) ______ +8.2 . 10-6 +0-6 . 10-6 tg.a . ^g-e 



error overestimates the uncertainty due to the truncation of the perturbative series. In the Hue of this 
discussion, and following [1], we will not use this prescription to estimate the systematic error on the 
truncation of the series, and we will limit the analysis to the uncertainties coming from the study of CIPT 
only. 

The discrepancies found between FOPT and CIPT at \s\ = are reduced drastically when computing 

i?^+^^_(M|) (see Fig. 6 and Table 2). The small value of as(M|) ensures a much better convergence 
of the perturbative series. The better convergence also leads to a tiny scale dependence, which is even 
smaller for CIPT than for FOPT, and hence to small theoretical uncertainties. 



3.3 Quark-Mass and Nonperturbative Contributions 

Following SVZ [38], the first contribution to Rr beyond the D = Q perturbative expansion is the non- 
dynamical quark-mass correction of dimension Z? = 2, i.e., corrections scaling like The leading 
D = 2 corrections induced by the light-quark masses are computed using the running quark masses 
evaluated at the two- loop level (denoted m in the following). The evaluation of the contour integral in 

FOPT [4] leads to terms oc m^^^(m^)/m^, TO„(m^)mrf(m^)/r7i^, which arc small. 

The dimension D = A operators have dynamical contributions from the gluon condensate {gsGG) and 
the light u,d quark condensates {niiqiqi), which are the vacuum expectation values of the gluon field 
strength-squared and of the scalar quark densities, respectively. The remaining D = A operators in- 
volve the running quark masses to the fourth power. Solving the contour integral [4] results in terms 

^^udv/A o.1i'm-T){0'sGG)/m^, {niqqq) /m^, 04(771^/771^), where remarkably the contribution from the 

gluon condensate vanishes at the first order in asim^). 

The contributions from dimension D ~ 6 operators are more delicate to analyse. The most important 
operators arise from four-quark terms of the form qj^riqjq^r2qi. We neglect other operators, such as the 
triple gluon condensate whose Wilson coefficient vanishes at order as, or those which are suppressed by 
powers of quark masses, in the evaluation of the contour integrals performed in [4]. The large number of 
independent operators of the four-quark type occurring in the D = 6 term can be reduced by means of 
the vacuum saturation assumption [38]. The operators are then expressed as products of (two-)quark con- 
densates as{^J-){qiqi{^J-)){qjqj{^J-))■ Since the scale dependence of the four-quark and two-quark operators 
are different, such factorisation can hold for a specific value of the renormalisation scale (at best). To take 
into account this problem as well as likely deviations from the vacuum saturation assumption, one can 
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introduce an effective parameter p (in principle scale-dependent) to replace the four-quark contribution 
by pasiqqy- The effective D = 6 term obtained in this way is [4] Sl^Jy^^^ oc pas{qq)'^ /m^, with a relative 
factor of —7/11 between vector and axial vector contributions. 

The D = 8 contribution has a structure of non-trivial quark-quark, quark-gluon and four-gluon conden- 
sates whose explicit form is given in [58] . For the theoretical prediction of Rr it is customary to absorb 
the whole long- and short-distance parts into the scale invariant phenomenological D — 8 operator (Os), 
which is fit simultaneously with as and the other unknown nonperturbative operators. Higher-order con- 
tributions from £> > 10 operators to Rr are expected to be small since, like in the case of the gluon 
condensate, constant terms and terms in leading order in as vanish after integrating over the contour. 
We will not consider these terms in the following. 

3.4 Impact of Quark-Hadron Duality Violation 

A matter of concern for the QCD analysis at the t mass scale is the reliability of the theoretical de- 
scription, i.e., the use of the OPE to organise the perturbative and nonperturbative expansions, and the 
control of unknown higher-order terms in these series. A reasonable stability test consists in varying ttIt 
continuously to lower values ^/sq < ttit for both theoretical prediction and measurement, which is possible 
since the shape of the full r spectral function is available. This test was successfully carried out [37, 25, 1] 
and confirmed the validity of the approach down to sq ^ 1 GeV^ with an accuracy of 1-2%. In this 
section, we consider a different test of the sensitivity of the analysis to possible OPE violations. 

The SVZ expansion provides a description of the correlator 77 (or of the Adler function D) for values 
of the incoming momentum in the deep Euclidean region, based on the separation between large and 
soft momenta flowing through the diagrams associated to this correlator. If the OPE description were 
accurate, we could check the cogency of this description by performing an analytic continuation of the 
OPE to any value of the momentum in the physical region and comparing it with the spectral functions in 
Fig. 2. As seen from these figures, perturbative QCD describes the asymptotic behaviour of the functions, 
but fails to reproduce their details. 

The OPE suffers from a similar failure as can be expected from the intrinsic nature of the OPE pro- 
cedure [38, 39, 40, 41, 42]: it only yields a truncated expansion in the first powers of l/Q, i.e., the 
singularities near x = of U'^'^ [cf. Eq. (3)). Therefore, it misses singularities for finite or — > oo 
related to long-distance effects. Even a large momentum q fiowing through the vacuum polarisation di- 
agrams may be split into a soft quark-antiquark pair and soft gluons: this physical possibility cannot 
be properly described by OPE, since no separation can be performed between hard and soft physics in 
such a situation. One expects for some of these effects to yield terms proportional to exp(— AQ)/(5'° or 
exp(— A^Q^)/Q^ (where are positive and A is a typical hadronic distance), which are exponentially 
suppressed in the deep EucHdean region and thus absent in the truncated OPE series. But once these 
terms are continued analytically along the branch cut, they generate a (power suppressed or exponentially 
suppressed) oscillatory behaviour of the spectral function, which is similar to the one in Fig. 2. Such a 
behaviour is generally called "violation of local quark-hadron duality" . 

To determine i?^., we compute the convolution of the OPE expression of the Adler function with a kernel 
along the circle of radius sq . We know that duality violation will have a small impact for the two regions 
close to the real axis (these terms are exponentially suppressed in the Euclidean region, and the kernel 
vanishes for s = sq). But to assess the systematic uncertainties related to the use of OPE, it is instructive 
— even if very approximate — to simulate the contributions of duality violating terms on the rest of 
the circle. For this purpose, we use two different models proposed in [41], which provide a coarse and 
rather qualitative description of such effects (one of these models has been very recently reconsidered 
in [44] to investigate duality-violating effects on the determination of nonperturbative condensates from 
ALEPH data in the vector channel). In both cases, one does not aim at a complete description of the 
correlator 77, but focuses on the deviation between the full description and its truncated OPE expansion 
Z\77 = 77 — 77*^^^. In the first model (7) the quarks propagate in an instanton background field with a 
fixed size p, leading to the duality violation 

An^'\Q) - ^K,iQp)K^,{Qp) , (27) 

where the -f'^(-)i are modified Bessel functions of the second kind. The second model (77) mimics a comb 
of resonances with a width that grows with the energy, so that they overlap progressively when the energy 



15 



increases 

"""'«> = -iTd?75^(*'W 4)- 

Here ip{z) is the di-gamma function, and z = {Q"^ / a'^Y~^/'^^ , where a parametrises the offset between 
the resonances, and B their (growing) widths. In this model, one can define 77'-'^^ as the expansion in 
powers of 1/z (up to 2;^ here, since we neglect operators of Z? = 10 and beyond). Duality violations are 
encoded in AU'^^^^ = Cii{n — n'^^^)'-^^K The factors C/j/ are normalisation constants. 

One can check that the two models share the same features: they are exponentially suppressed in the 
Euclidean region, and exhibit a branch cut for time-like values of s, such that they contribute to the 
spectral functions with oscillations decreasing in amplitude when the energy increases. They differ by 
the dependence of their oscillation frequency on the energy: the instanton model oscillates like sin(^/sp), 
while the resonance model varies like sin(s/(T). 

To investigate the numerical impact of quark-duality violation on our results, we vary for each model 
the parameters and fix the normalisation such that the imaginary part of sum of the perturbative QCD 
computation and of the duality-violating terms match smoothly the V + A spectral function near s = 
m^. We then compute the contribution of the duality- violating part to S^^^ by performing the contour 
integral (18). For the instanton model we asymptotically reproduce the data for p values between 2.4 and 
4.4 GeV~^, leading to a contribution to 5^"^ below 4.5 • 10"'^. For the resonance model we find values for 
between 1.65 and 2 GeV^, and B between 0.3 and 0.6, leading to a contribution to (5^"-' below 7 ■ 10^^. 
These limits are however quite conservative because the models used exhibit significant oscillations in 
the y -f A spectral function. Although allowed by the ALEPH data because of the larger error bars close 
to the endpoint, such oscillations are disfavored by the overall pattern of the spectral function, with 
oscillation amplitudes that are strongly suppressed above 1 GeV. Even though these two models could be 
improved in many ways, it is hard to see how their contributions to could be enhanced by an order of 
magnitude such that they would invalidate the OPE approach. At least in the case of the V + A spectral 
function, we therefore expect the violation of quark-hadron duality to have a negligible impact on our 
results. In the next section, we will see that the induced error on 

5(0) 

remains well within the systematic 

uncertainties coming from other sources. 



4 Combined Fit 

Apart from the perturbative term, the full OPE contains contributions of nonperturbative nature parametrised 
by higher-dimensional operators, whose value cannot be computed from first principles. It was shown 
in [5] that one can exploit the shape of the spectral functions via weighted integrals to obtain additional 
constraints on cks(m^) and — more importantly — on the nonperturbative power terms. 



4.1 Spectral Moments 

The r spectral moments at sq = "m-T defined by 



R^%A- I ds(l^^](^]:^^, (29) 



■.VIA - J \^^^) ) 



where R'^'^y/j^ — Rr.v/A- Using the same argument of analyticity as for R^-, one can reexpress (29) as a 
contour integral along the circle \s\ = so- The factor (1 — s/m'^)'' suppresses the integrand at s = 
where the validity of the OPE is less certain and the experimental accuracy is statistically limited. Its 
counterpart {s/rn^Y projects upon higher energies. The spectral information is used to fit simultaneously 
as(m^) and the leading D — 4,6,8 nonperturbative contributions. Due to the intrinsic experimental 
correlations (all spectral moments rely on the same spectral function) only four moments are used as 
input to the fit. 

In analogy to R^ (16), the contributions to the moments originating from perturbative QCD and non- 
perturbative OPE terms are separated. The prediction of the perturbative contribution takes the form 

00 

^(o.fe^) (30) 
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with the functions [1] 



1 

27ri 



ds 
s 



2r(3 + fc) 



r(i + £) 



r{2 + £) 



r{A + k + I) r{5 + k + £) 



J( — ,1 + ^,3 + fc 
So 



2/1 —,2 + £,3 + k 



(31) 



which make use of the elementary integrals I{'y,a,b) — Jq t"" ^(1 — t)^ ^dt. The contour integrals are 
numerically solved for the running as{—^s) using the CIPT prescription. 

In the chiral limit and neglecting the small logarithmic s dependence of the Wilson coefficients, the dimen- 
sion D nonperturbative contributions to the spectral moments simplify greatly (c/. matrix (133) 
in [1]). One finds that with increasing weight £ the contributions from low dimensional operators vanish. 
For example, the only nonperturbative contribution to the moment R].\^y^ 
and beyond (neglected). 



■"^^ stems from dimension D 



For practical purpose it is more convenient to define moments that are normalised to the corresponding 
Rt,v/a to decouple the normalisation from the shape of the t spectral functions. 



T^ke 



V/A 



(32) 



The two sets of experimentally almost uncorrelated observables — Rt.v/a oh one hand, and the moments 
-^tV/a ^^"^ other hand — yield independent constraints on Q!s(m^) and thus provide an important 
test of consistency. The correlation between these observables is negligible in the V + A case where 
Rry+A is calculated from the difference Rr — i?T,s, which is independent of the hadronic invariant mass 
spectrum. One experimentally obtains the D'^^y^j^ by integrating weighted normalised invariant mass- 
squared spectra. The corresponding theoretical predictions are easily adapted. 

The measured A and {V -\- A) spectral moments and their linear correlations matrices are given in 
Tables 3 and 4, respectively. Also shown are the central values of the theory prediction after fit convergence 
(c/. Sec. 4.2). The correlations between the moments are computed analytically from the contraction of the 
derivatives of two involved moments with the covariance matrices of the respective normalised invariant 
mass-squared spectra. In all cases, the negative sign for the correlations between the k = l,£ — and 
the k = 1,£ > 1 moments is due to the p (V) and the tt, ai (A) peaks, which determine the major 
part of the k = 1,£ — moments. They are less prominent for higher moments and consequently the 
amount of negative correlation increases with £ = 1,2,3. This also explains the large and increasing 
positive correlations between the fc = 1,^ > 1 moments, in which, with growing i, the high energy tail is 
emphasised more than the low energy peaks. The total errors for the {V + A) case are dominated by the 
uncertainties on the hadronic branching fractions. 



4.2 Fit Results 

Along the hue of the previous analyses from ALEPH [23, 37, 59, 1], CLEO [24], and OPAL [25], we 
simultaneously determine asirn^), the gluon condensate, and the effective D = 6,8 nonperturbative 
operators from a combined fit to i?r and the spectral moments D^^yj^^ with fc = 1, £ = 0, 1, 2, 3, taking 
into account the strong experimental and theoretical correlations between them. 

The fit minimises the of the differences between measured and predicted quantities contracted with the 
inverse of the sum of the experimental and theoretical covariance matrices. The theoretical uncertainties 
include separate variations of the unknown higher-order coefficient K^, for which the value/error if 5 = 
K4^{K/^/ K^) K, 378 ± 378 has been used, and of the renormalisation scale. The latter quantity has been 
varied within the range ± 2 GeV (corresponding to ^ = 1 ± 0.63), and the maximum variations of 
the observables found within this interval are assigned as systematic uncertainties (c/. Sec. 3.2.3). To 
avoid double counting of errors the estimated term has been fixed when varying ^. The corresponding 
systematic errors for as(m^) are 0.0062 {K^) and io^oo40 (0- The errors induced by the uncertainties on 
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Table 3. Experimental {DI.^y/a) ^-nd theoretical (D^ y^^° , obtained after fit convergence, cf. Sec. 4.2) spectral 
moments of inclusive vector {V), axial- vector (A) and vector plus axial- vector {V + A) hadronic r decays. The 
errors A'^^'^ D^y/A summarise statistical and systematic uncertainties. 







I = 1 


1 = 2 


^ = 3 


j-^li (theo) 
^exp^n^ 


0.71668 
0.71568 
0.00250 


0.16930 
0.16971 
0.00043 


0.05317 
0.05327 
0.00054 


0.02254 
0.02265 
0.00041 


^T,A 

j-^ie (thoo) 

T,A 

^cxp^lf^ 


0.71011 
0.71660 
0.00182 


0.14903 
0.14571 
0.00063 


0.06586 
0.06574 
0.00036 


0.03183 
0.03130 
0.00025 


^T,V+A 
p.ie (theo) 


0.71348 
0.71668 
0.00159 


0.15942 
0.15767 
0.00037 


0.05936 
0.05926 
0.00033 


0.02707 
0.02681 
0.00025 



Table 4. Experimental correlations between the moments v/a/v+a- Correlations between Rr,v+A, determined 
from the leptonic r branching fractions, and the corresponding moments are negligible. 

nio nil n^o n^^ nis n^^ n^s 

J^T,V ^T,V J-^T,V ^T,V ^T,A J^T,A ^T,A J^T,A T ,V + A J^T,V+A ^T,V+A 

Rr,v -0.287 0.153 0.274 0.302 -0.255 0.013 0.178 0.272 DI°v+a -0.722 -0.974 -0.987 

Dl°y 1 -0.821 -0.981 -0.993 Df^A 1 -0.746 -0.963 -0.978 D^^+a 1 0.801 0.662 

Dl]y - 1 0.899 0.824 D\]a " 1 0.866 0.646 D^^^+a " 1 0.975 
- - 1 0.988 D\^^A - - 1 0.938 



S'ew and \Vud\ amount to 0.0007 and 0.0005, respectively. With these inputs, the massless perturbative 
contribution (5*^*^^ is fully defined, and the parameter as(m^) can be determined by the fit. 

Table 5 summarises the results for the V , A and V A- A combined fits using CIPT. The (5(2) term is not 
determined by the fit, but is fixed from a theoretical input on the light quark masses varied within their 
errors [1]. The quark condensates in the S^'^^ term are obtained from partial conservation of the axial- 
vector current (PCAC), while the gluon condensate is determined by the fit, as are the higher-dimensional 
operators {Oq) and (Os). 

The advantage of separating the vector and axial-vector channels and comparing to the inclusive V + A 
fit becomes obvious in the adjustment of the leading nonperturbative contributions of Z3 = 6 and D = 
8, which have different signs for V and A and are thus suppressed in the inclusive sum. The total 
nonperturbative contribution, (5np — S'^'^'' +6^^'^ +S^^\ from the V+A fit, although non-zero, is significantly 
smaller than the corresponding values from the V and A fits, hence increasing the confidence in the ag (m^) 
determination from inclusive V + A observables. 

There is a remarkable agreement within statistical errors between the asirn^) determinations using the 

vector and axial-vector data, with Q;s^''(m^) — ai^^m^) — 0.013 ± 0.013, where the error takes into 
account the anticorrelation in the experimental separation of the V and A modes. This result provides 
an important consistency check since the two corresponding spectral functions are experimentally almost 
independent, they manifest a quite different resonant behaviour, and their fits yield relatively large 
nonperturbative contributions compared to the V + A case. Contrary to the vector case, the axial-vector 
fit has a poor value originating from a discrepancy between data and theory for the ^ = 0, 1 normalised 
moments (cf. Table 3). Although the origin of this discrepancy is unclear, it may indicate a shortcoming 
of the OPE in form of noticeable inclusive duality violation in this channel. The observed systematic 
effect on the asirn^) determination in this mode appears however to be within errors. From the fit to the 
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Table 5. Fit results for Qs(m^) and the nonperturbative contributions for vector, axial- vector and V+A combined 
fits using the corresponding experimental hadronic widths and spectral moments as input parameters, and using 
the CIPT prescription for the perturbative prediction. Where two errors are given the first is experimental and 
the second theoretical. The j'^-* term comes from theoretical input on the light quark masses varied within their 
allowed ranges (see text). The quark condensates in the 5'''' term are obtained from PC AC, while the gluon 
condensate is determined by the fit. The total nonperturbative contribution is the sum (5np = 

5(4) +5(6)^^(8)^ 



Parameter 


Vector (V) 


Axial- Vector (A) 


V + A 




0.3474 ± 0.0074t^ E!!j^^ 


0.3345 ± o.oo78i;;;!5J!?^ 


0.3440 ± 0.0046155;;;!!^;^ 


5(0) 

5(2) 


0.2093 ± 0.0080 
(-3.2 ±3.0) ■ 10"'' 


0.1988 ±0.0087 
(-5.1 ±3.0) • 10-'' 


0.2066 ± 0.0070 
(-4.3 ±2.0) ■ 10-* 


(a.GG) (GeV^) 

5(4) 
5(6) 

5(8) 

Total (5np 


(-0.8 ±0.4) ■ 10^^ 
(0.1 ±1.5)- 10-* 
(2.68 ±0.20) • 10-2 
(-8.0 ±0.5) ■ 10-=* 
(1.89 ±0.25) ■ 10-2 


(-2.2 ±0.4) • 10-2 
(-5.9 ±0.1) • 10-^ 

(-3.46 ±0.21) ■ 10-2 
(9.5 ±0.5) • 10-^ 

(-3.11 ±0.16) ■ 10-2 


(-1.5 ±0.3) ■ 10-2 
(-3.0 ±0.1) • 10-3 
(-3.7 ± 1.7) • 10-3 
(8.1 ±3.6) • 10-* 
(-5.9 ± 1.4) ■ 10-3 


xVdf 


0.07 


3.57 


0.90 



V + A T spectral function, we obtain 

as(m2) = 0.344 ± 0.005 ±0.007 , (33) 

where the two errors are experimental and theoretical. The values of the gluon condensate obtained in 
the V, A, and V + A fits are not very stable. Despite the apparent significance of the result for V + A, 
we prefer to enlarge the error taking into account the discrepancies between the V/A results. We find for 
the combined value {usGG) = (—1.5 ±0.8) • 10^^ GeV'', which is at variance with the usual values quoted 
in the applications of SVZ sum rules. We note however that not much is known from theoretical grounds 
about the value of the gluon condensate [57] . 

The result (33) can be compared with the recent determination [6], as(m^) — 0.332 ± 0.005 ± 0.015, also 
at N^LO, but using as experimental input only Rt^v+a, and not including the new information given in 
Sec. 2. Another major difference with our analysis is that both perturbative procedures, FOPT and CIPT, 
are considered on equal footing, and their results are averaged. This leads to the lower value for as{m'^) 
and to an inflated theoretical error including half of the discrepancy between the two prescriptions. 

The evolution of the value (33) to Aff , using Runge-Kutta integration of the four-loop /3-function [51], 
and using three-loop quark-flavour matching [62, 64, 65, 66], gives 

ai''^ {Ml) = 0.1212 ± 0.0005 ± 0.0008 ± 0.0005 , 

= 0.1212 ±0.0011 . (34) 

The first two errors in the upper line are propagated from the Q;s(m^) determination, and the last error 
summarises uncertainties in the evolution.^" All errors have been added in quadrature for the second line. 
The result (34) is a determination of the strong coupling at the Z-mass scale with a precision of 0.9%, 
unattained by any other as{Ml) measurement. The evolution path of as{m'^) is shown in the upper plot 
of Fig. 7 (the two discontinuities are due to the chosen quark- flavour matching scale of /i = Iniq). The 
evolution is compared in this plot with other as determinations compiled in [60] (we also included [63]), 
and with new NNLO measurements based on hadronic event shapes from e^e~ annihilation covering the 
energy range between 91.2 and 206 GeV [61]. 

*" The evolution error [1] receives contributions from the uncertainties in the c-quark mass (0.00020, rUc varied 
by ±0.1 GeV) and the 6-quark mass (0.00005, m;, varied by ±0.1 GeV), the matching scale (0.00023, fj. varied 
between Q.TWiq and 3.0mg), the three-loop truncation in the matching expansion (0.00026) and the four-loop 
truncation in the RGE equation (0.00031), where we used for the last two errors the size of the highest known 
perturbative term as systematic uncertainty. These errors have been added in quadrature. 
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1 10 10^ 

|i scale (GeV) 



Fig. 7. Top: The evolution of as(m^) to higher scales /i using the four-loop RGE and the three-loop matching 
conditions applied at the heavy quark-pair thresholds (hence the discontinuities at 2mc and 2rnb). The evolution 
is compared with independent measurements (taken from the compilation [60], and including the recent measure- 
ments [61, 63]) covering fi scales that vary over more than two orders magnitude. Bottom: The corresponding as 
values evolved to Mz- The shaded band displays the r decay result within errors. 



The theoretically most robust precision determination of as stems from the global fit to electroweak 
data at the Z-mass scale. As for as{m^), this determination benefits from the computation of the N'^LO 
coefficient occurring in the radiator functions that predict the vector and axial- vector hadronic widths 
of the Z (and also in the prediction of the total W width). We use the newly developed G/itter package [67] 
for the fit, and obtain 

af' {Ml) = 0.1191 ± 0.0027 ± 0.0001 . (35) 

The value and first error represents the fit result, and the second error is due to the truncation of the 
perturbative series. It is estimated similarly to the t case by adding a fifth-order term proportional to 
K^, estimated by K4^{K4^/ K:i)^ to the massless part, and a fourth-order term (estimated accordingly), 
containing large logarithms \a-{rrit / M z) : to the massive part. We also vary the renormalisation scale of 
the massless contribution within the interval ^ = 1 ± 0.63, assuming the fifth order coefficient to be 
known. The result (35) agrees with the finding of Ref. [6] . 

The r-based result (34) appears now twice more accurate than the determination from the Z width. Yet 
the errors are very different in nature with a t value dominated by theoretical uncertainties, whereas the 
determination at the Z resonance, benefiting from the much larger energy scale and the correspondingly 
small uncertainties from the truncated perturbative expansion, is limited by the experimental precision of 

the electroweak observables. The consistency between the two results, a's\^'z) ~ <^s^\Mz) ~ 0.0021 ± 
0.0029, provides the most powerful present test of the evolution of the strong interaction coupling as 
it is predicted by the nonabelian nature of QCD over a range of s spanning more than three orders of 
magnitude. The q^s\m1) determination agrees with the average of the three currently most precise full 
NN(N)LO measurements (deep inelastic scattering [68, 60], ALEPH event shapes between 91 and 206 
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GeV [61], and global electroweak fit at Mz), yielding an average of 0.1189 ± 0.0015 (0.1204 ± 0.0009) 
when not including (including) the r result, which is justifiably assuming uncorrelated errors. The r-based 
result differs at the 2.5cr level from the value 0.1170±0.0012 found in lattice QCD calculations with input 
from the mass splitting of the T resonances [69]. The average of all five values reduces the discrepancy 
to 2.1a (x^ probability of 0.04). 



5 Conclusions 

We have revisited the determination of q;s(to^) from the ALEPH r spectral functions using recently 
available results. On the experimental side, new BABAR measurements of the e+e^ annihilation cross 
section into KKir using the radiative return method now permit, through CVC, a much more accurate 
determination of the vector/axial-vector fractions in the corresponding r decays. Also, better results 
are available on t decays into strange final states from BABAR and Belle. On the theory side, the first 
unknown term in the perturbative expansion of the Adler function, the fourth-order term K4, was recently 
calculated, opening the possibility to further push the accuracy of the theoretical analysis of the hadronic 
r decay rate. 

Motivated by these improvements we have reexamined the theoretical framework of the analysis. In 
particular the convergence properties of the perturbative expansions for the t and Z hadronic widths 
have been studied, and the ambiguity between the fixed-order (FOPT) and contour-improved (CIPT) 
approaches for summing up the series has been discussed. The study confirms our earlier findings (at 
third order) that CIPT is the more reliable treatment. Furthermore we have identified specific consistency 
problems of FOPT, which do not exist in CIPT. Possible violations of quark-hadron duality at the r mass 
scale have been considered using specific models, and their effect has been found to be well within our 
quoted overall theoretical uncertainty (however, due to the coarseness of the models, we do not introduce 
additional theoretical errors). 

We perform a combined fit of the r hadronic width and hadronic spectral moments resulting in the value 
as(m^) = 0.344±0.005cxp±0.007thoo, consistent with the previous value obtained for three known orders, 
and with a 20% reduced theoretical uncertainty. This somewhat moderate improvement is the result of 
the relatively large value K4 ~ 49, suggesting a slowly converging perturbative series and giving rise to 
relatively large truncation uncertainties. Nevertheless, the result confirms the excellent accuracy that can 
be obtained from the analysis of r decays, albeit indicating that this method may approach its ultimate 
accuracy. 

The evolved r result at the Mz scale, Q;s(M|) = 0.1212 ± 0.0005oxp± O.OOOSthm ± O.OOOScvoi, is the most 
accurate determination available. It agrees with the corresponding value directly obtained from Z decays, 
which we have reevaluated. Both determinations are so far the only results obtained at N"^LO order. They 
confirm the running of as between 1.8 and 91 GeV as predicted by QCD with an unprecedented precision 
of 2.4%. 

We are indebted to Martin Gobel and the Gfitter group for implementing the new N^LO term into the global 
electroweak fit, and AH acknowledges the fruitful collaboration. We thank Oscar Cata, Maarten Golterman and 
Santi Peris for letting us preview their analysis on duality violation in hadronic r decays (which arrived after 
finalising this paper), and for helpful discussions. Many thanks to Matthias Jamin for pointing out a mistake in 
the total nonperturbative contributions previously quoted in Table 5 (corrected in the present version) . This work 
was supported in part by the EU Contract No. MRTN-CT-2006-035482, "FLAVIAnet". 
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